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^ ' Abstract 

^h ' Burgers turbulence subject to a force f(x,t) = J2j fj( x )d(t ~ tj)i where the tj's 

are "kicking times" and the "impulses" fj(x) have arbitrary space dependence, 



combines features of the purely decaying and the continuously forced cases. With 
large-scale forcing this "kicked" Burgers turbulence presents many of the regimes 
proposed by E, Khanin, Mazel and Sinai (1997) for the case of random white-in- 
time forcing. It is also amenable to efficient numerical simulations in the inviscid 
limit, using a modification of the Fast Legendre Transform method developed for 



c 

i-rt . decaying Burgers turbulence by Noullez and Vergassola (1994). For the kicked 

case, concepts such as "minimizers" and "main shock" , which play crucial roles in 
recent developments for forced Burgers turbulence, become elementary since ev- 
erything can be constructed from simple two-dimensional area-preserving Euler- 
Lagrange maps. 

The main results are for the case of identical deterministic kicks which are 

/\^ periodic and analytic in space and are applied periodically in time. When the 

space integrals of the initial velocity and of the impulses vanish, it is proved and 
illustrated numerically that a space- and time-periodic solution is achieved ex- 
ponentially fast. In this regime, probabilities can be defined by averaging over 
space and time periods. The probability densities of large negative velocity gra- 
dients and of (not-too-large) negative velocity increments follow the power law 
with —7/2 exponent proposed by E et al. (1997) in the inviscid limit, whose ex- 
istence is still controversial in the case of white-in-time forcing. This power law, 
which is seen very clearly in the numerical simulations, is the signature of nascent 
shocks (preshocks) and holds only when at least one new shock is born between 
successive kicks. 

It is shown that the third-order structure function over a spatial separation Ax 
is analytic in Ax although the velocity field is generally only piecewise analytic 
(i.e. between shocks). Structure functions of order p / 3 are nonanalytic at 
Ax = 0. For even p there is a leading-order term proportional to |Ax| and for 
odd p > 3 the leading-order term oc Ax has a nonanalytic correction oc Ax|Ax| 
stemming from shock mergers. 



1 Introduction 

The driven Burgers equation 

d t u + ud x u = vd 2 x u + /, (1) 

u(x,t ) =u (x), (2) 

with a force f(x,t) has a much richer structure than the decaying problem with / = 0. 
Indeed, by the Hopf (1950) and Cole (1951) transformation the latter is mapped into the 
heat equation, while the former goes into a kind of imaginary-time Schrodinger equation 
with a potential F such that / = —d x F. Actually, the randomly forced Burgers equation 
appears in a number of problems in statistical mechanics. This includes the Kardar, Parisi 
& Zhang (1986) equation for interface dynamics (see also Barabasi & Stanley 1995) and the 
problem of directed polymers in random media (Bouchaud, Mezard & Parisi 1995). 

From the point of view of turbulence dynamics, the forced Burgers equation can be used 
in the same spirit as the forced Navier-Stokes equation, namely to investigate universality 
with respect to the forcing of various statistical properties. For Navier-Stokes turbulence, 
when the force is confined to large spatial scales and the Reynolds number is very high, 
small-scale (inertial range) statistical properties are generally conjectured not to depend on 
the forcing, except through overall numerical factors. Similar conjectures have been made for 
Burgers turbulence with large-scale forcing. For example, there is little doubt that, because 
of the presence of shocks, all the structure functions of order p > 1 have a universal exponent 
equal to unity (see, e.g., Bouchaud, Mezard & Parisi 1995; E et al. 1997). Chekhlov & Yakhot 
(1995) were the first to study the behavior of the probability density function (pdf) of velocity 
increments when the force is a white-noise process in time and has a k _1 spectrum, a case 
leading to anomalous scaling. When the force is a white-noise process in time but is confined 
to large scales, the behavior of the pdf p(£) of the velocity gradient £ at large negative values 
is rather controversial. If it is generally believed that this pdf follows a power law, 

p(0<x\$\ a , for£^-oo, (3) 

the conjectured values of a differ markedly. Polyakov (1995) and Boldyrev (1997), using a 
field-theoretical operator product expansion, predict a = —5/2; E et al. (1997), using a semi- 
heuristic approach in which nascent shocks (preshocks) are key, predict a = —7/2; Gotoh & 
Kraichnan (1998), using a Fokker-Planck equation approach, predict a = —3; more recent 
work by Kraichnan (1999) favors a = —7/2. E & Vanden Eijnden (1999, 2000) develop a 
probabilistic formalism that copes with the delicate problems arising in the limit of vanishing 
viscosity when shocks are present, prove that a < —3 and make a good case for a = —7/2. 
The question of the correct law for the case of white-noise forcing remains however open (we 
shall come back to this in §||). There are simpler situations for which the arguments in favor 
of a = —7/2, originally developed by E et al. (1997), can be made rigorous, such as one- and 
multi-dimensional decaying Burgers turbulence with smooth random initial conditions (Bee 
& Frisch 2000; Bee, Frisch & Villone 2000). 

It may be thought that numerical experimentation on the one-dimensional forced Burg- 
ers equation should be able to easily obtain the correct scaling laws. This is actually rather 
difficult if one tries to use standard numerical schemes of the kind also applicable to the 
Navier-Stokes equation, such as spectral methods with dissipation explicitly taken into ac- 
count. Indeed, it has been shown by Gotoh & Kraichnan (1998; see also Gotoh 1999) that, 
in the presence of a small but finite viscosity, there is a range of large negative values of the 
gradient for which the pdf p(£) oc ^|£| _1 , which decreases rather slowly at large |£|. This 



is a direct consequence of the hyperbolic-tangent internal structure of shocks. As a result, 
the behavior of the pdf at negative values of £ smaller in absolute value is contaminated and 
may display power laws somewhat shallower than predicted by the theory in the inviscid 
limit. See, for example, figure 3 of Gotoh & Kraichnan (1998) in which a power-law range 
with exponent -3 is seen over a little more than one decade of £- values. Such artefacts will 
disappear if much higher resolution is used (Gotoh and Kraichnan were using between 2 17 
and 2 20 collocation points). It is, however, much more efficient to use alternative numerical 
schemes where one works directly with the inviscid limit. We shall see that such schemes, 
which were previously used only for the decaying case can be readily extended to the forced 
case and are particularly well suited for the case of kicked Burgers turbulence, in which the 
force is concentrated at discrete times. 

The paper is organized as follows. Section || introduces the general concept of kicked 
Burgers turbulence and gives an explicit representation of the solution which is a simple 



generalization of the so-called "minimum representation" for the decaying case. Section 2.1 



presents the concepts of "minimizers" and "main shocks" for the case when the dynamics are 



started at to = ~~ °°- Section |2.2j presents the Fast Legendre Transform numerical scheme for 
space-periodic kicked Burgers turbulence with spatially smooth forcing. The next two sections 
are about the case of deterministically kicked Burgers turbulence when the kicks are periodic 
in space and time. Section || presents numerical results on exponential convergence to a 
periodic solution and shows how this is related to properties of minimizers (the rigorous results 
on convergence to a unique periodic solution are derived in the Appendix). Section |4| presents 



the main results about pdf's of velocity derivatives and increments (§4.1) and about structure 
functions (§f^). Section g presents concluding remarks and some possible extensions to 
Navier-Stokes turbulence. 

2 The inviscid limit for kicked Burgers turbulence 

We shall be concerned here with the initial-value problem for the one-dimensional Burgers 
equation (|l|) when the force is concentrated at discrete times: 

/M) = £/j(x) *(*-*;), (4) 

3 

where 5 is the Dirac distribution and where both the "impulses" fj(x) and the "kicking 
times" tj are prescribed (deterministic or random). The kicking times are ordered and form 
a finite or infinite sequence. In this paper the impulses are always taken smooth and acting 
only at large scales. Newman & McKane (1997) have used similar kicking, but confined to 
small scales, in a context where the forced Burgers equation is used for the study of directed 
polymers. Kraichnan (1999) has considered a simple model in which there are non-smooth 
impulses creating directly sawtooth profiles in the velocity. The precise meaning we ascribe 
to the Burgers equation with this forcing is that at time tj, the solution u(x,t) changes 
discontinuously by the amount fj{x) 

u(x,tj+) = u(x,tjJ) + fj(x), (5) 

while, between tj+ and t(j+i)- the solution evolves according to the unforced Burgers equa- 
tion. Without loss of generality, we can assume that the earliest kicking time is tj = to, 
provided we set fj = uo and u(x, t) = for t < to- 

It is clear that any force f(x,t) which is continuously acting in time can be approximated 
in such a way by choosing the kicking times sufficiently close. 



We shall also make use of the formulation in terms of the velocity potential ip(x,t) and 
the force potentials Fj{x) 



u(x,t) = -d x ip(x,t), fj(x) = -—Fj(x). 



The velocity potential satisfies 



dtip = - (d x i(j) 2 + vd xx ^ + Y, Fj(x) S(t - tj), 

3 

ip(x,t ) = ipo(x), 



(6) 

(7) 
(8) 



where ipo(x) is the initial potential. 

As it is well-known, the solution to the unforced Burgers equation with positive viscosity 
is, has an explicit integral representation obtained by Hopf (1950) and Cole (1951) which 
can be used to investigate the limit of vanishing viscosity. We are here exclusively inter- 
ested in this limit. Generically, shocks appear then after a finite time from smooth initial 
data. (The correct solution may also be obtained by solving the inviscid equation with a 
variational formulation and the condition that, at a shock, the right velocity is less than the 
left velocity (Lax 1957; Oleinik 1957).) Use of Laplace's method then leads to the following 
"minimum representation" for the potential in the limit of vanishing viscosity (henceforth 
always understood) which relates the solutions at any two times t > t' between which no 
force is applied: 

(x - yf 



ip(x,t) 



mm 

y 



2(t - t> 



i>(y,t') 



0) 



It is known that, when if is the initial time, the position y which minimizes @ is the 
Lagrangian coordinate associated to the Eulerian coordinate x. The map y i— ► x is called the 
Lagrangian map. By expanding the quadratic term it is easily shown that the calculation 
of ip(-,t) from ifj(-,t') is equivalent to a Legendre transformation. For details on all these 
matters, see She, Aurell & Frisch (1992) and Vergassola et al. (1994). 

We now turn to the forced case with impulses applied at the kicking times tj. Let tju\ 
be the last such time before t. (Henceforth we shall often just write tj.) Using @ iteratively 
between kicks and changing the potential ip(y,tj+i) discontinuously by the amount Fj + \(y) 
at times tj+i, we obtain 



ip(x,t) 

Mh;x,t;{yj}) 



min [A(jo;x,t]{yj})-i/>o(y jo )], 



(a - yj) 2 

2(t - tj) 



j-i 

+ E 

3=30 



2(tj+i - tj) 



F j+i(Vj+i) 



(10) 

(11) 



where A(j ; x,t; {yj}) is called the action. 

For the Burgers equation with a continuous-in-time force deriving from a potential F(x, t), 
E et al (1997, 2000) give a minimizer representation of the solution: 



ip(x,t) 



A(t ;x,t;y(-)) = 



min [A(t ;x, t; y{-)) - ipo(y(to))] 



to 



f{s) 



F(y(s),s) 



ds. 



(12) 
(13) 



where the minimum is taken over all curves y(s) satisfying y(t) = x. Note that this repre- 
sentation, which actually goes back to work by Oleinik (1957) on general conservation laws, 
is just a continuous limit of (|Io|) -(PT|), obtained by taking tj = j At and letting At — > 0. 



Returning to the case of kicked Burgers turbulence, from (|l0|), we shall now introduce 
the concept of "minimizers" and "main shock" (§ |2.1| ). Eq. ( JlOJ ) will also be our starting 
point for the numerical method (§ |2.2|) . For the rest of the paper, we shall assume that the 
force potential and the initial condition are periodic in the space variable. For convenience, 
the period is taken to be unity in the theory, while 2-7r-periodicity is assumed in numerical 
studies. 

2.1 Minimizers and main shocks 

For the case of the kicked Burgers equation with an initial condition at tj a "minimizing 
sequence" associated to (x, t) is defined as a sequence of yj's (j = jo, jq + 1,..., J(t)) at which 



the r.h.s. of ( JlQj ) achieves its minimum. Differentiating the action (11) with respect to the 
yj's one gets necessary conditions for such a sequence, which can be written as a sequence of 
(Euler-Lagrange) maps 

Vj+i = Vj + fj(yj), (14) 

Vj+i = Vj + v j+ i(t j+1 -tj) = yj + (vj + fj(Vj))(t j+ i -tj), (15) 



where 



Vj ~ Vj-i 



■> ~ i — i 

bj bj-\ 

These equations must be supplemented by the initial and final conditions: 



(16) 



v h = Mv j ), (17) 

x = yj + v J+ i(t-tj). (18) 

It is easily seen that u(x,t) = vj + \ = (x — y,j)/{t — tj). Observe that the "particle velocity" 
Vj is the velocity of the fluid particle which arrives at yj at time tj and which, of course, has 
remained unchanged since the last kick (in Lagrangian coordinates). Eq. ( |i~4] ) just expresses 
that the particle velocity changes by fj(yj) at the the kicking time tj. 

Note that (14)-(|l5[) define an area-preserving and (explicitly) invertible map. 



The presence of a force, deterministic or random, allows a formulation of the Burgers 
equation in the semi-infinite time interval ] — oo, t] without fully specifying the initial condition 
uq{x) but only its (spatial) mean value (u) = J uo(x)dx. Heuristically, this follows from the 
observation that, for a force of zero spatial mean value, as assumed here, (u) is a first integral, 
and hence does not depend on time, while all the other information contained in the initial 
condition is eventually forgotten. 

Actually, the construction of the solution in a semi-infinite time interval is done by ex- 
tending the concept of minimizing sequence to the case of dynamics starting at to = — oo. 



For a semi-infinite sequence {yj} (j < J), let us define the action A(— oo; x, t; {yj}) by (11) 
with jo = — oo. Such a semi-infinite sequence will be called a "minimizer" (or "one-sided 
minimizer") if it minimizes this action with respect to any modification of a finite number of 
yj's. Specifically, for any other sequence {yj} which coincides with {yj} except for finitely 
many j's (i.e. yj = yj, j < J — k, k > 0), we require 

A(J - k; x, f, {yj}) > A(J - k; x, t; {%•}). (19) 

Of course, the Euler-Lagrange relations (14)-(15) still apply to such minimizers. Hence, if 



for a given x and t we know u(x,t) we can recursively construct the minimizer {yj} backwards 
in time by using the inverse of (|l4])-(15) for all j < J and the final condition - now an initial 



condition - (|lg|) with v j + \ = u{x, t). This is well defined except where u(x, t) has a shock and 
thus more than one value. Actually, solutions in a semi-infinite time interval are constructed 
from minimizers and not the other way round. 

One way to construct minimizers is to take a sequence of initial conditions at different 
times to ~~ ► — °°- At each such time some initial condition uq(x) is given with the only con- 
straint that it have the same prescribed value for (u). Then, (finite) minimizing sequences 
extending from to to t are constructed for these different initial conditions. This sequence 
of minimizing sequences has limiting points (sequences themselves) which are precisely min- 
imizers (E et al. 2000). The uniqueness of such minimizers, which would then imply the 
uniqueness of a solution to Burgers equation in the time interval ] — oo, t], can only be shown 
by using additional assumptions, for example for the case of time-periodic forcing (§H and 
Appendix). 

If (u) = 0, the sequence {yj} minimizes the action A{— oo;x, t; {yj}) in a stronger sense. 
Consider any sequence {yj} such that, for some integer P we have jjj = yj+P, j < J—k, k > 
and which differs arbitrarily from {yj} for j > J — k. (In other words, in a sufficiently remote 
past the hatted sequence is just shifted by some integer multiple of the spatial period.) We 
then have 

A(-oo;x,t;{yj}) > A(-oo;x, t; { Vj }). (20) 

Indeed, for (u) = 0, the velocity potential for any initial condition is itself periodic. In this 
case a particle can be considered as moving on the circle S 1 and its trajectory is a curve 
on the space-time cylinder. The yj's are now defined modulo 1 and can be coded on a 
representative < yj < 1. The Euler-Lagrange map (|l4|)-(^) is still valid provided (|l5|) is 
defined modulo 1. 

The condition of minimality implies now that yj and yj+i are connected by the shortest 
possible straight segment. It follows that \vj + \\ = p{yj,yj + i)/{tj + \ — tj), where p is the 
distance on the circle between the points yj,yj + i, namely p(a,b) = min{|a — b\,l — \a — b\}. 
Hence, the action A can be rewritten in terms of cyclic variables: 



A(-oo;x,t;{yj}) 



P 2 (x,yj) 
2(t - tj) 



+ £ 



j<JL 



p 2 (yj+i,yj) 

2(t j+1 - tj) 



Fj+iiv. 



j+i/ 



(21) 



We now introduce the concept of "global minimizers" (or "two-sided minimizers" ) limiting 
ourselves to the case (u) = for simplicity. We first observe that any minimizer {yj,j < 
J} can be continued for all j > J and hence times t' > t by using the system (|i~4|)-(|l5|). 
However, this procedure, when extended too far in time, will not usually generate a minimizer 
associated to time t'. Nevertheless, for any time t there always exist positions x such that the 
corresponding minimizers {yj, j < J(t)} can be continued to the bilateral sequence {yj, — oo < 
j < +00} while keeping the minimizing property. Such global minimizers correspond to 
trajectories of fluid particles which, from t = —00 to t = +00, have never been absorbed in 
a shock. 

We then observe that any shock existing at time t can be continued for all times s > t : 
shocks can merge but they cannot otherwise disappear. However, since new shocks can be 
produced, it is not always possible to trace back an existing shock for arbitrary times s < t. 
A shock with this property of having always existed in the past is called a "main shock" . 

In E et al. (2000) it is shown that for the case of random forcing which is 1-periodic 
in space and white noise in time rather than impulsive, the solution of Burgers equation in 
] — 00, t] is unique for (u) = 0. It is also shown that at time t, the set of points x G S 1 
with more than one minimizer, that is shock locations, is finite and that the main shock 
and the global minimizer are unique. The global minimizer forms a hyperbolic trajectory of 



the Euler-Lagrange equations and all other minimizers approach the global one, as t — > — oo 
exponentially fast. In particular the two minimizers associated to the main shock approach 
the global minimizer in the remote past but it may be shown that they do so from opposite 
directions on the circle S . We shall see below in §|3| that the same picture holds in the case 
of generic time-periodic kicking. 

2.2 A Fast Legendre Transform numerical method 

The numerical method used here solves the kicked Burgers equation directly in the inviscid 
limit. The basic ideas are very simple: at each kicking time tj the potential is changed by 
the amount Fj(x); between two successive kicks (or between the last kick and the output 
time) the decaying Burgers equation is solved using the minimum representation (||); this 
procedure is repeated as many times as the number of kicks between the initial time and the 
output time. 

Specifically, the space periodic interval, here taken to be [0, 2n[, is discretized on a regular 
grid of N collocation points x^ = 2kTt/N. For each of these positions, (||) is used to determine 
the potential at a time t in terms of the potential just before the last kicking time tj, time 
at which the potential has discontinuously changed by Fj{x). Hence, we have : 



i>k(t) = tf>{x k ,t) = -min 



%^-Wt,J+F,( w) ) 



(22) 



where ipt(tj-) is the value of the velocity potential at x = ye, just before the kick. We note 
yi k the minimizing position corresponding to Xk- This procedure is, in principle, applied 
recursively, starting from tj . 

The problem is that naive application of (|22j ) yields an algorithm with 0(N 2 ) operations 
between two successive kicks. She, Aurell &: Frisch (1992) observed that the minimizing 
position y is actually a monotonic non-decreasing function of x. This is indeed a simple con- 
sequence of the convexity of the parabolic term involved in (^). Hence, the determination of 
the ye^s can be performed using a binary-subdivision search which requires only 0{N dog 2 N) 
operations. This kind of algorithm is known under the name Fast Legendre Transform (FLT), 
since the minimum representation is equivalent to a Legendre transform (Hopf 1950). (An 
even faster algorithm requiring only 0{N) operations has been developed by Trussov (1996).) 

We use an adaptation of the method of Noullez and Vergassola (1994; see also Vergassola 
et al. 1994) who developed an FLT algorithm using a binary-subdivision search combined 
with a reorganization of the search, permitting the use of very low in-core storage. We first 
determine the minimizing yi Q for the point xq = 0. As the velocity-potential at the time tj + is 
periodic, it is easy to show that ye is within the interval [— tt, +vt[. The minimizing location 
corresponding to xn = 2ir is then given by periodicity and reads ye N = yi + 2ir = yi +N- 
The search for all the other minimizing locations yi k can then be restricted to indices i^ such 
that £o < Ik < £q + N. We then compute yi N/2 , corresponding to x = it. We can then further 
subdivide the x-interval by considering k = N/4 and k = 3N/4, for which the corresponding 
4's satisfy £ < ?n/4 < ^-N/2 < 4iV/4 < 4 + N. At the next stage, we compute yi N/s , yi 3N/8 , 
yi 5N/s and yi 7N/s for which we need £q, £n/4, £n/2 an d £zn/a as search boundaries. We repeat 
this subdivision procedure log 2 N times to obtain the iV values of yi k . 

The method just described is optimal for non-smooth solutions of the kind considered 
by Vergassola et al. (1994) who had initial conditions of Brownian type. For the case 
of smooth solutions considered here, a more accurate determination of the solution is re- 
quired to obtain reliable results on space derivatives of the velocity. We now describe an 
improvement of the method allowing to calculate first- and second-order derivatives. We 



observe that in (|22|), when the discrete location y% is replaced by an arbitrary real num- 
ber y, the minimum, for a given xj~ is, in general, not achieved exactly on the grid at 
yi k , but at a neighboring location y{k) within less than one mesh. This location satisfies 
%k = y(k) + (t — tj) [u(y(k),tj-) + fj(y(k))], obtained by requiring that the derivative of 
the r.h.s. of ( |22|) vanish. For this, the velocity u(y,tj-) and the force fj(y(k)) are Taylor 
expanded to second order to obtain the improved location y{k). When required, the first and 
second space derivatives of the velocity u(x, t) are calculated, not from the potential by finite 
differences, but by using exact expressions of these derivatives in terms of the Lagrangian 
map from the preceding kicking time, which holds for the unforced Burgers equation (E et 
al. 1997; Bee & Frisch 1999). 

This FLT, implemented on a grid of 2 17 ~ 10 5 collocation points takes about 1 s of CPU 
on a 100 MFlops computer. Without use of the binary-subdivision search the CPU time 
would be several thousands times larger. Among the other advantages of FLT is that no 
viscosity is needed and that the solution can be calculated directly at the required output 
times without need to obtain it at many intermediate times (other than the kicking times). 

3 Deterministic periodic kicking 

From now on we shall consider exclusively the case where the kicking is periodic in both 
space and time. Specifically, we assume that the force in the Burgers equation is given by 

+oo 

f(x,t) = g(x) J2 ^-JT), (23) 

j=-oo 

g(x) = -^G(x), (24) 

where G(x), the kicking potential, is a deterministic function of x which is periodic and 
sufficiently smooth (e.g. analytic) and where T is the kicking period. The initial potential 
V'init^) is also assumed smooth and periodic. This implies that the initial velocity integrates 
to zero over the period. (The case where this assumption is relaxed will be considered briefly 
in the Conclusion in relation with the Aubry-Mather theory.) 

The numerical experiments reported hereafter have been made with the kicking potential 

G{x) = - sin3x + cosx, (25) 

and a kicking period T = 1. Other experiments where done with (i) G(x) = —cosx and 
(ii) G(x) = (l/2)cos(2x) — cosx. The former potential produces a single shock and no 
preshock. As a consequence it displays no —7/2 law in the pdf of gradients. (In this paper, 
we limit ourselves to cases presenting at least one new shock between successive kicks.) The 
latter potential gives essentially the same results as reported hereafter but has an additional 
symmetry which we avoided by the choice (p5|). 

The number of collocation points chosen for our simulations is generally N x = 2 ~ 
1.31 x 10 5 , with a few simulations done at N x = 2 20 (for the study of the relaxation to the 
periodic regime presented below) . Since our numerical method allows us to go directly to the 
desired output time (from the nearest kicking time) there is no need to specify a numerical 
time step. However, in order to perform temporal averages, e.g. when calculating pdf's or 
structure functions, without missing the most relevant events (which can be sharply localized 
in time) we need sufficiently frequent temporal sampling. We have taken for the total number 
of output times N t ~ 1000 chosen such that the increment between successive output times 




Figure 1: Snapshots of the velocity for the unique time-periodic solution corresponding to 
the kicking force g{x) shown in the upper inset; the various graphs correspond to six output 
times equally spaced during one period. The origin of time is taken at a kick. Notice that 
during each period, two new shocks are born and two mergers occur. 



is roughly the two-thirds power of the mesh (this is related to the structure of preshocks, see 



0). 



Figure [I] shows snapshots of the time-periodic solution at various instants. It is seen that 
shocks are always present (at least two) and that at each period two new shocks are born at 
ti,\ ~ 0.39 and t+2 ~ 0.67. There is one main shock which remains near x = tv and which 
collides with the newborn shocks at t c \ ~ 0.44 and t c i ~ 0.86. Figure |3| shows the evolution 
of the positions of shocks during one period. 

We find that, for all initial conditions uo(x) used, the solution u(x, t) relaxes exponentially 
in time to a unique function 1*00(2;, t) of period 1 in time. Figure ^ shows the variation of 
Jq 77 \u(x, n_) — u OQ (x, 1_)| dx/(2Tr) for three different initial conditions as a function of the 
discrete time n. 

The phenomenon of exponential convergence to a unique space- and time-periodic solution 
is something quite general: whenever the kicking potential G(x) is periodic and analytic and 
the initial velocity potential is periodic (so that the mean velocity (u) =0 at all times), there is 
exponential convergence to a unique piecewise analytic solution. This is proved rigorously in 
the Appendix for functions G{x) which have a unique point of maximum with a nonvanishing 
second derivative (Morse generic functions). Here, we just explain the main ideas of the proof 
and give some additional properties of the unique solution. 

One very elementary property of solutions is that, for any initial condition of zero mean 
value, the solution after at least one kick satisfies |tt(x,i)| < (l/2) + maXa; \dG(x)/dx\. Indeed, 
at a time t = n_ just before any kick we have x = y+u(x, n__) where y is the position just after 
the previous kick of the fluid particle which goes to x at time n_. It follows from the spatial 
periodicity of the velocity potential that the location y which minimizes the action is within 
less than half a period from x. Thus, \u(x, n_)| < 1/2. The additional max x \dG(x)/dx\ term 
comes from the maximum change in velocity from one kick. It follows that the solution is 
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Figure 2: Exponential relaxation to a time-periodic solution for three different initial velocity 
data as labelled. The horizontal axis gives the time elapsed since t = 0. 
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Figure 3: Evolution of shock positions during one period. The beginnings of lines correspond 
to births of shocks (preshocks) at times £*i and £*2! shock mergers take place at times t& and 
t C 2- The "main shock", which survives for all time, is shown with a thicker line. 
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Figure 4: Sketch of a hyperbolic fixed point P with stable (r' s **) and unstable (r^ u ^) manifolds. 
The dashed line gives the orbit of successive iterates of a point near the stable manifold. 

bounded. Note that if the spatial and temporal periods are L and T, respectively, the bound 
on the velocity becomes L/(2T) + max x \dG(x) / dx\. 

The convergence at large times to a unique solution is related to properties of the two- 
dimensional conservative (area-preserving) dynamical system defined by the Euler-Lagrange 
map ([H|)-(|l5D of §2T. By construction, we have u(x, 1+) = u(x) — dG(x)/dx, where u(x) is 



the solution of the unforced Burgers equation at time t = 1_ from the initial condition u(x) 
at time t = 0+. The map u h-» u(x) +g(x), where g(x) = —dG(x)/dx, will be denoted B g and 
is just the map which solves the kicked Burgers equation over a time of one. The problem is 
to show that the iterates B'^uq converge for n — ► oo to a unique solution. 

If it were not for the shocks it would suffice to consider the two-dimensional Euler- 
Lagrange map. Note that, for the case of periodic kicking, this map has an obvious fixed 
point P, namely (x = x c ,v = 0), where x c is the unique point maximizing the kicking 
potential. It is easily checked that this fixed point is an unstable (hyperbolic) saddle point 
of the Euler-Lagrange map with two eigenvalues A = 1 + c + \/c 2 + 2c and 1/A, where 
c=-(l/2)(d 2 G(x c )/dx 2 ). 

Like for any two-dimensional map with a hyperbolic fixed point, there are two curves 
globally invariant by the map which intersect at the fixed point: the stable manifold r( s ), the 
set of points which converge to the fixed point under indefinite iteration of the map, and the 
unstable manifold TW, the set of points which converge to the fixed point under indefinite 
iteration of the inverse map, as illustrated in figure |] (see, e.g., Henon 1983, Manneville 1990). 
It follows that any curve which intersects the stable manifold transversally (the tangents of the 
two curves are distinct) will, after repeated applications of the map, be pushed exponentially 
against the unstable manifold at a rate determined by the eigenvalue 1/A. In the language 
of Burgers dynamics, the curve in the (x,v) plane defined by an initial condition Uq(x) will 
be mapped after time n into a curve very close to the unstable manifold. In fact, for the 
case studied numerically, 1/A ~ 0.18 is within one percent of the value measured from the 
exponential part of the graph shown in figure ^. Note that if the initial condition Uq(x) 
contains the fixed point, the convergence rate becomes (1/A) (even higher powers of 1/A are 
possible if the initial condition is tangent to the unstable manifold) . 

The fixed point P gives rise to a very simple global minimizer: (jjj = x c , Vj = 0) for all 
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Figure 5: Unstable manifold r' u ' on the (x, u)-cylinder (the ^-coordinate is defined modulo 1) 
which passes through the fixed point P = [x c , 0). The bold line is the graph of u oc (x, 1_). 
The main shock is located at x\ = x r . Another shock at x\ corresponds to a local zig-zag of 
T (u) between A and B. 

positive and negative j's. It follows indeed by inspection of ([y) that any deviation from 
this minimizer can only increase the action; actually, it minimizes both the kinetic and the 
potential part of the action. Note that the corresponding fluid particle is at rest forever 
and will never be captured by a shock (it is actually the only particle with this property). 
It is shown in the Appendix that any minimizer is attracted exponentially to such a global 
minimizer as t — ► — oo. Thus, any point (jjj,Vj) on a minimizer belongs to the unstable 
manifold r^ n > and, hence, any regular part of the graph of the limiting solution Uoo(x) belongs 
to the unstable manifold r' u '. This unstable manifold is analytic in the relevant region but 
can be quite complex. It can have several branches for a given x (see figure |5|) and does not 
by itself define a single- valued function Uoo(x). The solution has shocks and is only piecewise 
analytic. Consideration of the minimizers is required to find the position of the shocks in the 
limiting solution: two points with the same x corresponding to a shock, such as A and B on 
figure [5] should have the same action. 

Finally, we give the geometric construction of the main shock, the only shock which 
exists for an infinite time. Since A is positive, locally, minimizers which start to the right of 
x c approach the global minimizer from the right, and those which start to the left approach 
it from the left. Take the rightmost and leftmost points x r and xi on the periodicity circle 
such that the corresponding minimizers approach the global minimizer from the right and left 
respectively (see figure |6|). These points are actually identical since there cannot be any gap 
between them that would have minimizers approaching the global minimizer neither from the 
right nor the left. The solution u^x) has then its main shock at xi = x r . 
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Figure 6: Minimizers (trajectories of fluid particles) on the (x, £)-cylinder. Time starts at 
— oo. Shock locations at t = 0_ are characterized by having two minimizers (an instance is 
at x\). The main shock is at xi = x r . The fat line x = x c is the global minimizer. 

4 Statistical properties for the periodically kicked case 

We are here working with time- and space-periodic deterministic solutions of the kicked 
Burgers equation. We thus choose to define our statistical averages as averages over the two 
periods, here assumed to be both unity. Specifically, let F{u) be an observable (functional 
of the solution u) and let T x t denote the space-time translation operator which shifts the 
solution u by a spatial amount x and a temporal amount t. We define 



(T(u))= f f T(T X:t u)dxdt. 

Jo Jo 



(26) 



For example, with the observable J-(u) = [u(Ai, 0) — u(0, 0)] p , we obtain the structure func- 
tion of order p over a separation Ax: 



S p (Ax)= I / [u(x + Ax,t) -u(x,t)] p dxdt. 
Jo Jo 



(27) 



Such averages are easily calculated numerically. For example, pdf's are obtained from 
space-time histograms over all collocation points and a suitably large number of output 
times. 

4.1 Pdf's of velocity derivatives and increments 

For the periodic solution of §H we calculate first and second space derivatives of the veloc- 
ity; the corresponding pdf's are then determined as normalized space-time histograms after 
binning of derivative values (the bins are in geometric progression ; there are 100 bins per 
decade for the first derivative and 50 for the second). Figures |7| and || show the pdf's of the 
first and second space derivatives in log- log coordinates. Negative values are shown for the 
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Figure 7: Pdf of the velocity gradient at negative values in log-log coordinates. Upper inset: 
local scaling exponent. A power law with exponent —7/2 is obtained at large arguments. 

former and positive values for the latter. It is seen that clean power laws are obtained. More 
quantitative information about the values of the exponents of the power laws are obtained by 
measuring the "local scaling exponent" , i.e. the logarithmic derivative of the pdf, calculated 
here using least-square fits on quarter decades. The results are shown as upper right insets 
on figures and g[ It is seen that, over one decade, the local exponent for the pdf of the 
gradient is within less than five percent of the value —7/2 predicted by a simple theoretical 
arguments given hereafter; for the second derivative there are about four decades within five 
percent of the value —2. 

The presence of a power-law pdf with —7/2 exponent is easily understood. It is just the 
signature of the preshocks which appear when new shocks are created during the unforced 
phase between two successive kicks. These preshocks are the only structures giving large 
finite negative gradients: shocks give infinite negative gradients (unless a finite viscosity 
is introduced) and the gradients in the immediate spatial neighborhood of a mature shock 
are not particularly large. The theory of the —7/2 law at large negative values £ of the 
velocity gradient, developed by Bee & Frisch (2000) for the unforced case with random initial 
conditions is readily adapted to the deterministic case, provided we use space-time averages 
instead of ensemble averages. A simplified presentation, following in part E et al. (1997), 
is given hereafter for the case of a single preshock. The contributions of several preshocks 
to the pdf are just additive and it is proved in the Appendix that the periodic solution has 
finitely many preshocks. 

We define a velocity in Lagrangian coordinates u^' (a, t') with the origin of time just after 
a kick and a = at a (negative) minimum of d x u. Without loss of generality, we assume 
vS '(0, 0) = (otherwise we perform a Galilean transformation to bring it to zero). We then 
have, locally, u( L )(a,0) = — c\a + C2<2 3 + h.o.t., where c\ and C2 are positive constants and 
"h.o.t." stands for higher-order terms. No generality is lost by assuming c\ = 1 (otherwise 
make a linear change on the a-coordinate). The fluid particle initially at a will be at time 
t' at x = a + t'vS '(a, 0) = a(l — t') + t'c2a? + h.o.t. This "Lagrangian map" becomes 
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Figure 8: Same as figure with the second space derivative of the velocity. The exponent is 
now —2. 

singular at t 1 = t* = 1, the instant of preshock (formation of a shock). We then have 
a = {x/c2) 1 ' i + h.o.t. (The cubic root is here defined both for positive and negative values 
of its argument.) Since the Lagrangian velocity has not changed, the Eulerian velocity is 
given by k(x,£*) = — (X/C2) 1 ' 3 + h.o.t., which has a cubic root structure and a gradient 
— (x/c2) _2 ' 3 /(3c2). Hence, the gradient takes large negative values for small x. Just before 
£*, at time t' = 1 — r, we have, x = ra + C2<i 3 + h.o.t. It follows that the cubic relation 
between a and x still holds, except in a region of Lagrangian width ~ r ' and thus of 
Eulerian width ~ r 3 ' 2 , where the relation becomes linear to leading order. (It is because of 
this t 3 ' 2 dependence that the time r between successive outputs and the mesh Sx = 2tt/N 
must be related by 5x ~ r 3 / 2 .) 

The question is now: what is the fraction of space-time where the velocity gradient 
d x u < £i where £ is a large negative number? Because of the cubic root structure, x must 
be in a small interval of width ~ |£|~ 3 • The time must be sufficiently close to £* for this 
interval still to be in the region of validity of the cubic relation, that is, within ~ |x| 2 ' 3 ~ |£| _1 - 
Hence, the relevant space-time fraction or, in other words, the cumulative probability to have 
d x u < £ is ~ |£l - • This gives a pdf ~ |£|~ 7 ' 2 at large negative £'s. 

Actually, there is another contribution, also proportional to |£|~ stemming from a small 
time interval r ~ |x| 2 ' 3 ~ |£| just after t* when small-amplitude shocks are present which 
have not yet completely destroyed the cubic root structure (Bee & Frisch 2000). For the case 
studied numerically, where the kicking potential is given by (^3), there are two preshocks, 
each giving a contribution to the pdf of the gradient proportional to |£|~ 7 ' 2 . 

This argument is readily adapted to second space derivatives, yielding a pdf ~ |£|~ 2 as 
observed in figure |8|. (The same law holds also at large positive values since the second 
derivative near a preshock is an even function.) 

We now turn to the pdf of (spatial) velocity increments over a separation Ax. We define 



Au(Ax; x, t) = u(x + Ax, t) — u(x, t). 



(28) 



Its pdf's for various values of Ax are again calculated from space-time histograms. One 



15 




Figure 9: Upper part: Pdf of (negative) velocity increments in log-log coordinates for various 
values of the separation Ax in geometric progression from 2ir/N to 2 7 (2ir/N). Lower part: 
the corresponding local scaling exponents. 

hundred bins per decade are used. Ax is given the values 2ir2 p /N, where 2tt/N is the 
numerical mesh and p is varied from zero to seven. Figure pi gives log-log plots of the pdf 's of 
increments for the eight separations chosen. We limit ourselves to negative increments. The 
corresponding local scaling exponents are shown in figure ||. It is seen that, for moderately 
large increments, the pdf's have the same power-law behavior with exponent —7/2 as the pdf 
of the gradient. This universal behavior was also predicted by E et al. (1997) for white-in- 
time forced Burgers turbulence. Phenomenologically, this range is obtained simply by Taylor 
expanding the increment as Ax d x u. At larger increments (in absolute value) the local scaling 
exponent rises quickly to positive values but does not saturate to the value +1 predicted by 
E et al. (1997) by the following argument, based on the consideration of nascent shocks and 
which applies also to the periodically kicked case: The probability (as fraction of space) 
to have a shock in an interval of length Ax is oc Ax. Since the shock amplitude grows as 
(t — i*) 1 ' 2 , where t+ is the time of the preshock, the fraction of time for which the shock 
amplitude does not exceed a value |Au| is proportional to (Au) 2 . Hence, the cumulative 
probability to have a velocity increment (in absolute value) less than \Au\ is proportional to 
Ax(Au) 2 and the pdf is proportional to Ax|A-u|. By equating the contributions from the 
—7/2 and the +1 ranges, the transition between the two ranges is predicted to happen around 
an increment Au c which scales as (Ax) 1 ' 3 , in good agreement with our data. A clean +1 
range is not seen and would require a resolution of well over one million collocation points. 
Finally, the flat range seen in figure |] for \Au\ <C Ax, is a universal contribution oc (Ax) -1 
from extremal points of the velocity (not predicted by E et al. (1997)). This range extends also 
to positive values of Au <C Ax but there is no other universal range for positive increments; 
this is why pdf's are not shown for such increments. 
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Figure 10: Structure functions S p (Ax) for p = 2,3, 4, 5 as labelled. Note the linear behavior 
at small Ax. 



4.2 Structure functions 

We now study the structure functions for the limiting (unique) solution Uoo(x, t) of the period- 
ically kicked Burgers equation. For numerical studies 27r-periodicity in space and 1-periodicity 
in time are assumed. Hence, the structure function of (integer) order p is given by 



Y rl pin 

S p (Ax) = — dt dx[u oa (x + Ax,t) - u 00 (x,t)] t 
2tt Jo Jo 



(29) 



The 27r-periodicity of Uoo(x,i) immediately implies that S p (Ax) is 27r-periodic in Ax and is 
an even/odd function for even/odd p. 



Figure 10 shows the structure functions of order 2, 3, 4 and 5 (as labelled) in linear 
coordinates. It is seen that all these structure functions behave proportionally to Ax at 
small arguments (more precisely as Ax (sign (Ax)) p ). This is a well-known consequence of 
the presence of shocks (E et al. 1997). In the next two sections we shall show that all the 
structure functions except S^Ax) are nonanalytic functions of Ax. 



4.2.1 Analyticity of the third-order structure function 

For notational convenience in this and the next sections we assume 1-periodicity in space 
and time. Space averages over the period are denoted {•} x . Averages over both space and 
time are denoted (•). We shall prove that, when the kicking potential G(x) is analytic, the 
third-order structure function is also analytic. 



17 



This will be established as a consequence of the following relation for the unforced Burgers 
equation with space-periodic solution in the limit of vanishing viscosity. Let u = u(x, t) and 
v! = u(x + Ax, t), we have 

dt(u'u) x = ±d Ax ((u'-uf) x . (30) 

It is here assumed that Ax is not an integer multiple of the spatial period and that no pair 
of shocks remains separated by exactly Ax for a finite amount of time (this holds for almost 
every Ax). 

Proof. Let us denote by Xi(t) (i = 1, . . . , N(t)) the (Eulerian) ordered positions of shocks 

and by \[u]\i = u(Xi(t)+,t) — u(Xi(t)-,t) the (negative) velocity jump at the i-th. shock. 

(N(t) may change in time.) Except at shocks we can use the inviscid Burgers equation 

dtu + ud x u = 0. At shocks this has to be supplemented by the Rankine-Hugoniot conditions 

(see, e.g. Lax 1957) 

u 2 
Xi\[u]\i-\[-]\ i = ^ (31) 

which follow also from momentum conservation applied to small intervals straddling the 
shocks. An immediate consequence is that the velocity of the shocks are given by 

X i (t) = ^[u(X i (t)+,t)+u(X i (t)-,t)}. (32) 

We observe that 

(u'u) =} / u'udx, (33) 

where Xjv+i(i) = X\(t) + 1. Differentiating with respect to t, we obtain 

N 

d t {u'u) x = {ud t u' + u'd t u) x -Y J Xi{t)\Hiu{X i (t) + Ax,t) (34) 

4 = 1 



N rX i+ i(t) N 

J2 / {uu'd x u + u'ud x u) dx -J2 X i(t)\[u]\i uiX.it) + Ax, t) (35) 

1 J J\4 it) 1 



n Xi+l{t )n ! x 

}^ / t;uoaxU - -u 6axU dx 



t=i 



N N 2 

-Y^M^HiuiXiit) + Ax,t) + J2l\]\iuiXiit) + Ax,t). (36) 

j=l i=\ 



In going from (34) to (|35|) we used the inviscid decaying Burgers equation; from (pq) to 



we have performed an integration by parts and used d x v! = dj\ x u' . From (plj) follows that 



the last two terms in (36) cancel. Hence, we obtain 



dt (u'u) x = ^8 Ax (-uu >2 + u 2 u') x = ^8 Ax ((«' - uf) x , (37) 

which completes the proof. 

We now return to the case of the periodically kicked Burgers equation, with the unique 
solution Uoo(x,t). Using (|30|), integrated in time between two successive kicks, say at t = 
and t = 1, we have 

1~ ~ ,. . 1- / l3 



-dAxSsiAx) = -dAx([uooix + Ax,t) -^t 00 (a;,t)] , " 
6 6 \ 



(u OQ ix + Ax,l_)u 00 (x,l_)) x . - (u 00 (x + Ax,0 + )u oc ix,0 + )) x . (38) 
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Next, we use 

u QO (x, 0+) = Uoo(x, 0_) + g(x) = Uoo(x, 1_) + g(x), (39) 



which follows from (0) and (p3|); here, g(x) = —dG(x)/dx where G(x) is the kicking potential. 
Substituting this in (p£J), we obtain 

-a A x^ 3 (Ax) = (g(x)g(x + Ax)^ - ( 5 (x) Uoo (x + Ax, 0+))^ - (g(x + Ax) Uoo {x, 0+)) x . (40) 

We now assume that the kicking potential and, hence, g{x) are analytic functions and we 
find that all three terms on the r.h.s. of ( f40[) are analytic functions of Ax. This follows 
indeed from the observation that the analyticity of g(x) and the boundedness of h(x) imply 
the analyticity in Ax of the integral J Q g{x)h{x + Ax) dx, which is basically a convolution 
integral. We have thus proved the analyticity of the third-order structure function in the 
separation Ax. 

When the kicking potential G(x) has only a finite number of Fourier harmonics a stronger 
result holds: the third-order structure function has exactly the same harmonics as the kicking 
potential. This follows because the r.h.s. of ( ^0[) is a convolution integral. For the case of the 
kicking potential given by (|25|), which has the harmonics of wavenumber 1 and 3, we thus 
have 

S 3 (Ax) = A sin(Ax) + /i sin(3Ax). (41) 

(The presence of only sine functions is due to the odd character of the third-order structure 
function.) We have indeed checked that the structure function S3 (Ax) calculated numerically 



at the beginning of § |4.2| has a global fit of this form with A ~ —10.9953 and /i ~ —1.1463 
with an error of less than 10~ 5 . 

We finally observe that the analyticity result for the third-order structure function is quite 
general and has been proved also for the case of white-noise forcing (E & Vanden Eijnden 
2000). 

4.2.2 Nonanalyticity of the structure functions of order p 7^ 3 

We now concentrate on integer values of p > 1. Indeed, for non-integer values, the structure 
function is not defined, unless we take the absolute value of the velocity increment which 
results trivially in nonanalyticity and, for p = 1, the structure function vanishes. We intend 
to show that 

<? rAiO = f A p\ Ax \ + B p( Ax ) 2 + ° (i Ax ) 2 ) ' for even p; , , 

pl J 1 A p Ax + B p Ax\Ax\ + ((Ax) 2 ) , for odd p, [ ' 

where the constant A p never vanishes and the constant B p vanishes for p = 3 and never 
vanishes for p > 3 (the expressions of these constants will be given below). This will then 
imply (i) that all structure functions are proportional to the first power of the separation 
(a well-known result; see, e.g. E et al. 1997) and (ii) that all structure functions of order 
p / 3 are nonanalytic functions of Ax. Actually, we shall establish (|42|) only for Ax > 0; 
the extension to Ax < follows then from the even/odd character of structure functions of 
even/odd orders. 

The idea of the proof is to observe that the only possible sources of nonanalyticity are 
singularities of the solution in the space-time domain, namely, preshocks, shocks and shock 
mergers. The contributions from the analytic regions to S p (Ax) is clearly O ((Ax) p ) and must 
therefore be retained in (p2| ) only for p = 2. Let us now concentrate on the contributions 
from singularities. 
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It is easily shown that preshocks contribute at most terms O f (Ax)( p+5 '' 3 j which are 
only higher order corrections to (|42|). This follows from the scaling properties of the pdf of 
increments as discussed in § |4.1| and in Bee & Frisch (1999) which is itself a consequence of 
the cubic root structure of preshocks. 

As to the contribution of (mature) shocks, we obtain it by first calculating the contribu- 
tion, denoted St (Ax), coming from the neighborhood of individual shocks, ignoring shock 
mergers; then we determine the correction due to mergers, denoted AS*™ erge (Aa;). 

Let Xj(t) denote the positions of the various shocks (their number may change in time). 
Let iu(x, i) and uj(x,t) denote the velocity in the immediate right and left neighborhood 
of the j-th shock. Let Cj(t) = u^(Xj(t),t) — uj(Xj(t),t) denote the (negative) jump at 
Xj(t). (The jump Cj(t) is taken equal to zero when a shock has not yet been born or has 
disappeared by merger.) Since we take Ax > the requirement that x and x + Ax straddle 
the j-th shock limits the domain of x-integration to the interval ]Xj(t) — Ax,Xj(t)[. Hence, 
we have 



s; hock {Ax) 



dt 



j 



Xj{t) 
Xj(t)-Ax 



dx 



u 



~(x + Ax,t) 



U J 



(x,t) 



~\p 



(43) 



Since u^(x, t) and u- (x, t) are smooth functions of x, we can Taylor expand them near Xj(t). 
For our interest only the first two terms are relevant. We thus obtain 



S; hock {Ax)= dtJ2 



Xj(t) 



dx\uUX j (t),t) + d x uj (Xj(t),t) (x + Ax-Xj(t)) 



uJ (Xj(t),t) - \d x uj (X,(t),t)} (x - Xjit)) + o {(Ax) 2 ) . (44) 



We then use the following relation which governs the evolution of shock jumps: 

Cj(t) 



d_ 
dt 



T,Cj(t) 



d x uUx j (t),t) + d x uj(x j (t),t) 



(45) 



(This relation is obtained by using the inviscid Burgers equation on both sides of the shock 
and (£3^)-) Using ( |45|) in ([44|), performing all the space integrals and keeping only terms up 
to O ((Ax) 2 ), we obtain 



^ hock (Ax) =J o dtJ2 [c](t)Ax - -^-1 (Cf-\tj) (Ax) 2 ] + o ((Ax) 2 ) 



(46) 



We turn to the contributions of the finite set {-Mfc} of shock mergers taking place in the 
periodic space-time domain. (We only consider mergers of two shocks since events with more 
than two shocks merging are not generic; furthermore, it may be checked that they do not 
change our conclusions.) Associated with each event M&, we define: t^, the time of merger, 
Xki(t) and Xkr(t), the positions of the left and the right shocks about to merge and Cki(t) 
and Ckr(t), the respective jumps across these merging shocks. Let 9k be the first instant of 
time t when the distance Xk r (t) — X^i(t) becomes less than Ax. It is clear that, for small 
Ax, we have ik — 6k = O(Ax) (see figure |llj) . Furthermore, we have 



X kr (t) - X M (t) =Ax+(t- 9k) C ^k) + C kr (t k ) + o{Axy 



(47) 



After 6k, an interval ]x, x+Ax[ may either straddle a single of the two merging shocks or both. 
The calculation above did not take into account the possibility of straddling two shocks. (For 
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space 



Figure 11: Merger of two shocks Xki{t) and X\. r (t) on the (x, £)-cylinder. The merger (event 
Mfc) takes place at t = i&. At the earlier time t = Ok the shocks are within a distance Ax. 



Ok < t < tk, this happens when Xk r (t) — Ax < x < Xki(t).) It is also necessary to modify the 
spatial integration domain associated to situations where the interval of length Ax straddles 
a single shock. All these effects together give a correction to (|4q) stemming from mergers 
which, to leading order, reads 

(Cki(i k ) + Ckr{i k )Y - C p ki(ik) ~ Cl r {tk) 



A5 mcrge (Ax) = _^ x 



k 



Ckl{tk) + Ckr(tk) 



+ 



o((Ax) 2 ). (48) 



We now assemble the various contributions. The second term in the integral on the r.h.s. 
of (|46| ) can be integrated explicitly as a sum of terms coming from the birth and death of 
shocks. Preshocks do not contribute because they have vanishing jumps. Shock mergers 
give three contributions: two from the incoming shocks with jumps Cki{tk) and Ckr{tk) and 
one from the merged shock with jump C^(tfe) + Cf. r (tf.). We finally obtain the following 
expressions for the coefficients in the expansion ( fi"2| ) of the structure functions in terms of 
the shock jumps Cj(t): 
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(Ckl(tk) + Ckr(tk)) 
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C ll^k) 



(Ckl(tk) + Ckr(tk)) 



Cl(ik) 



Ckl(tk) + Ck r (tk) 



for all p > 2. 



(50) 



For p = 2, we must add to the expression given by the r.h.s. of (pOQ the contribution from 
the analytic regions which give also terms O ((Ax) 2 ). It is readily seen that B p given by ( |50| ) 
vanishes for p = 3 and only for that value. 

It must be stressed that, for structure functions of order p > 3, the oc (Ax) 2 corrections 
to the leading oc Ax terms come entirely from mergers. For cases which have a single shock 
with no mergers and no preshocks (an instance is G{x) = — cosx) this correction is absent. 



5 Concluding remarks 

We have shown that Burgers turbulence in the inviscid limit with periodic large-scale kicking 
is characterized by universal properties originally conjectured to hold for the case of random 
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forcing. In particular, there is a power-law tail with exponent —7/2 in the pdf of negative 
velocity gradients. This law was proposed by E et al. (1997) for the case of random forcing 
which is smooth in space and white noise in time. The validity of the —7/2 law for the latter 
case is still an open question. It is clear, that the —7/2 law is unescapable as soon as preshocks 
are present and well separated. Some hypothetical clustering of preshocks could invalidate 
the —7/2 law for the white- in-time case. Careful numerical experimentation using a sequence 
of random forces which approach white-in-time forcing should be able to shed some light 
on this issue. This can in principle be done using kicked Burgers turbulence in which the 
instants of kicking are taken closer and closer and the successive spatial impulses are taken 
random and independent. Obtaining sufficiently clean power-law scaling to distinguish, e.g. 
between an exponent —3 and —7/2 in the random case may require very large computational 
resources. It is also useful to investigate the statistical distribution of preshocks to test their 
possible clustering properties. 

Let us briefly now address the question of the effect of a finite small viscosity. Basically, 
this will broaden the shocks giving them a hyperbolic tangent structure of width oc v. From 
this it is easily inferred that the maximum negative gradient is |£| m ax = 0(V -1 ) and that the 
shoulders of such viscous shocks contribute a term oc ^|i;| _1 to the pdf of (negative) gradients 
(Gotoh & Kraichnan 1998). This term will dominate over the inviscid contribution oc |£|~ 7 ' 2 
beyond a crossover value of the gradient |£| c oc v~ 2 ' 5 . A small viscosity will also regularize 
preshocks, giving them a finite velocity gradient |£| m axpr oc u^ 1 ' 2 (Crighton &; Scott 1979). 
This gives only subdominant contributions for all £'s. 

The Burgers equation constitutes a dissipative dynamical system because of the presence 
of shocks which introduce an essentially irreversible element into the dynamics. Many features 
of our periodically kicked Burgers problem are actually in exact correspondence with those of 
a conservative (Hamiltonian) dynamical system, namely the equilibrium positions of a one- 
dimensional chain of (classical) atoms connected by elastic springs in the presence of a space- 
periodic external potential (Frenkel-Kantorova model). This problem has been investigated 
by Aubry (1983) and Mather (1982). The potential energy which has to be minimized to 
obtain the ground state has the following form: 

#({%}) = E \(Vi+i ~ Vi ~ a ? ~ eG (*)> ( 51 ) 

3 

where the yj's are the positions of the atoms, a is the unstretched length of the springs and 
e > measures the depth of the periodic external potential eG(x). Nontrivial properties of 
the ground states reflect the competition between the tendency of the atoms to sit at the 
minimum of the potential — eG{x) and to be within a distance a from each other. It is easily 
checked that the action (in the sense of §0) for the kicked Burgers equation with periodic 



forcing is exactly given by the Aubry-Mather Hamiltonian (|5lj) if we take a forcing potential 
eG{x) and a mean velocity (u) = a. The velocity in the Burgers equation is now the analog 
of the distance between adjacent atoms. Note that j is a space index in the chain model, 
whereas it is a time index in the Burgers equation. In this paper we have assumed a vanishing 
mean velocity (a = 0). In this case there is no competition between the aforementioned two 
tendencies, and the global minimizer y? = x c corresponds to a trivial ground state minimizing 
both parts of the Hamiltonian. More delicate effects arising in the case (u) ^ and in more 
than one dimension will be discussed in forthcoming work. We note that the connection 
between Aubry-Mather theory and Burgers equation for the case of time-periodic potentials 
with continuous time was discussed for the first time in Jauslin, Kreiss & Moser (1997) and in 
E et al. (2000). The theory has been recently developed further in E (1999) and in Sobolevski 
(1999). 
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Finally, we discuss some possible extensions of the ideas of the present paper to Navier- 
Stokes turbulence. Obviously, the method of forcing by kicks applied at discrete instants 
of time can also be used for Navier-Stokes. One of the things which made this method 
particularly valuable for the Burgers case is the existence of an efficient numerical algorithm 
to solve the purely decaying Burgers equation in the inviscid limit. At the moment there is 
nothing comparable for Navier-Stokes. Note that periodic kicking will not result in a unique 
time-periodic solution for Navier-Stokes since periodicity-breaking bifurcations leading to 
chaos will unavoidably occur. 

The results concerning analyticity of structure functions are likely to be the same for 
Burgers turbulence and isotropic 3-D Navier-Stokes turbulence (in the limit of vanishing vis- 
cosity). Indeed, in the latter case all structure functions of order p ^ 3 are generally believed 
to have scaling properties at small separations with nontrivial (and certainly noninteger) 
exponents. As to the third-order structure function, our proof of analyticity for analytic- 
in-space forcing can be extended if one assumes that Kolmogorov's four-fifths law is valid 
(Kolmogorov 1941; see also Chapter 6 of Frisch 1995). 

Burgers turbulence presents an algebraic tail for the pdf of the velocity gradient but 
nothing similar is known for Navier-Stokes turbulence. This may be telling us something 
about the possible singularities associated to the Navier-Stokes and Euler equations. A 
remarkable feature for the Burgers equation is that only preshocks but not mature shocks 
contribute to the power-law tail. There is indeed a basic difference between the two types 
of singularities. For analytic forcing, near a mature shock, the solution is not only piecewise 
analytic (i.e. on each side of the shock) but uniformly so: the radius of convergence of the 
Taylor series remains finite as one approaches a mature shock. In contrast, when approaching 
a cubic root preshock singularity, the radius of convergence goes to zero and gradients become 
very large. It is the algebraic behavior of their strength which causes the power-law tail. If 
the Navier-Stokes equation, in the limit of vanishing viscosity, were to develop any singularity 
of this kind (accompanied by algebraically large gradients) it should also display a power- 
law-tail pdf. Further increases in the quality of experimental and numerical turbulence or 
convection data are needed to find if such singularities are really ruled out. 
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APPENDIX 

Statement of the results for periodic kicking 

Here, we formulate and prove formally the statements presented somewhat heuristically in 

§y. The kicking is assumed 1-periodic in both space and time. The force in the Burgers 

equation is given by 

d + °° 

f(x,t) = -—G(x) J2 W-j), (A.l) 

j=-oo 

where G{x) is a deterministic function of x which is 1-periodic and three times continuously 
differentiable (G 6 C 3 ). For some of the statements below it is assumed that G is analytic. 
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We shall also assume that the kicking potential G(x) is generic in the Morse sense. This 
implies that G(x), considered on the circle < x < 1, has a unique point of maximum x c 
and that G(x) is non-degenerate at x c , i.e. d 2 G{x c )/dx 2 < 0. Without loss of generality we 
can assume that G(x c ) = 0. We denote by c = — (l/2)d 2 G(x c )/dx 2 > 0. The initial potential 
V'o(x) is also assumed 1-periodic. This implies that J uo(x)dx = 0, where uq{x) = u(x,0 + ) 
is the initial velocity. 

We now solve the unforced Burgers equation between the times t = and t = 1 and get 
u(x) = u(x,l-). Then, we determine u(x, 1+) = u(x) + g(x), where g{x) = —dG{x)/dx. 
Denote by B g the transformation from u(x,0+) to u(x, 1+) : B g u = u + g. Clearly, B g 
transforms Lq = {u(x) G i 1 [0, 1], J Q u(x)dx = 0} into itself. Then, the following statements 
hold. 

51 The functional transformation B g has a unique fixed point u^: B g Uoo = Uqo, 
max x |ttoo(x)| < 1/2 + max x |5(x)|. The fixed point Uqq is a function of bounded vari- 
ation, it is continuous everywhere except at the set of shock points, which is at most 
countably infinite. 

52 Let sbea point of continuity for Uoq. Then, for all uq such that J uo(x)dx = 0, 

BgUo(x) —> u 00 (x) as n — > oo. (A. 2) 

53 The unique global minimizer 7 C corresponds to a particle with zero velocity sitting at 
the point x c of maximum kicking potential. 

54 There exists a unique entropy weak solutionQn 00 (x, t) to the kicked Burgers equation in 
the semi- infinite domain ] — oo,T] with zero mean velocity. This solution is 1-periodic 
in time and it is generated by u^x): u OQ (x,k+) = Uoo(x) for all integer k. The 
solution Uoo(x,t) satisfies the estimate: max Ii( |uoo(x,t)| < 1/2 + max x |g(x)|. 

55 For arbitrary t there exists a unique main shock. 



S6 Convergence in ( |A.2| ) is exponentially fast in n. If Uoo is differentiable at x, then there 
exists a constant C(x) which does not depend on uq such that 

|£>oO) - uoo(x)| < C(x)\- n , (A.3) 



where A = 1 + c + \/c 2 + 2c > 1. If x c is a point of continuity for uq and uq{x c ) ^ 
then there exists a constant c{uq) > such that for all x the following estimate holds: 

\B n g u Q (x) - Uoo (x)\ > c(u )\- n . (A.4) 

S7 If the kicking potential G(x) is analytic, then for any t the solution Uoo(x, t) is a piece- 
wise analytic function of x. The number of pieces is finite and is equal to the total 
number of shocks at time t. The number of preshock events between time and 1 is 
also finite. 



1 Entropy solutions are the limit of viscous solutions as v — > (Lax 1957). 
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Proof 

Uniqueness of the solution 

Denote by T the Euler-Lagrange diffeomorphism of the phase-space cylinder C = {0 < y < 



1, — oo < v < 00} generated by the system of equations (15), i.e. T(y,v) = (y',v'), where 



v' = v + g(y) (A.5) 

y' = y + v + g(y) (modi). 



Since G G C 3 , the diffeomorphism T £ C 2 . The global minimizer in the sense of § |2.l| is a 
sequence {yj = x c } for all j. Here, it is trivial and it corresponds to a fixed point P = (x c , 0) 
for the diffeomorphism T, i.e. the corresponding trajectory of T is a stationary trajectory 
(yj,Vj) = (x c ,0) = P. Let {y'j, j < J, y'j = x} be an arbitrary minimizer. Then y' — ► x c 
as j — > —00. Indeed, if it were not the case, then a sequence {g/j = x, y" = x c , jq < j < 
J— I, y" = y'j, j < jo — 1} would have smaller action than {y'A for any jo which is sufficiently 
negative. (The action is given by (f2l|).) Convergence y'j — ► x c implies that the corresponding 
trajectory of T, {(y',f'- = ±p(y',y'_ 1 ))} converges to F as j — ► —cxd. Here, p denotes the 
distance between two points on the circle. Easy calculation shows that P is a saddle point 
for the diffeomorphism T with the eigenvalues Ai = A = l + c + \Jc 2 + 2c > 1 and < A2 = 
1/A < 1. It follows that there exist two C 2 -smooth curves r' s ) and T^ u ' which are stable and 
unstable manifolds for the point P. Both curves pass through P and consist of all points 
(y,v) whose trajectories approach P as j tends to +00 and —cxd, respectively (see figure |j). 
In fact, convergence is exponentially fast and the rate is given by 1/A. Since (y'-,t>') — > P 
as j — > —cxd, a point (y'j, v'j) belongs to r' u ^. Let us associate a minimizing curve 7 Xi j on 
the space-time cylinder M = {0 < y < 1, —cxd < t < +00} with an arbitrary minimizer 
{y'j, j < J}- To construct j x> j we just connect all the pairs of points (y'j,j), (y'j_\,j — 1) for 
j < J by straight segments of minimal length. Denote by y' x j(t) a piecewise linear function 
such that j Xt j = {(y' x j(t),t), —cxd < t < +cxd}. Clearly, y' x j(j) = y'j. Denote also by 7 the 
minimizing curve {(x c ,t), —00 <t< +cxd} corresponding to a global minimizer {yj = x c }. 

We have shown above that any minimizing curve y x j is asymptotic to 7 as t — ► —00. By 
a standard argument this implies that any two minimizing curves do not intersect each other, 
except if they start from the same point x. It follows that for all but at most count ably many 
x and all J there exists a unique minimizing curve j x ,j- Denote by S the exceptional set 
of x's where a minimizer is not unique. Obviously, for fixed x and different J's minimizing 
curves 7^7 are connected by a time shift. It follows that the set S does not depend on 
J. Hence, for all x outside of S one can define a function v(x) = dy' x j(t)/dt \t=j_ = 
±p(y'j,y' J _ 1 )- It is easy to see that for any solution u(x,t), t < T to the kicked Burgers 
equation in the semi-infinite domain ] — 00, T] one has: u(x, &_) = v(x) for all integer k < 
T and x outside of S. This implies that the solution to the kicked Burgers equation in 
the semi- infinite domain ] — 00, £] is unique and is generated by v(x). It also follows that 
"Uoo^) = v(x) + g(x) is a unique fixed point for B g . The set S is a set of shocks at integer 
moments of time. It follows from the closeness of the set of minimizers that v(x) is continuous 
outside of S. The non-intersecting property implies that for arbitrary x G S there exist 
v(x-) = lim y ^ x _v(y), v(x + ) = lim y ^ x+ v(y), and that v(x^) > v(x + ). It also follows 
that v (x) is a function of bounded variation. Clearly, v (x_) and v(x + ) are the velocities 
of two minimizing curves which start at (x, J) for any integer J. Notice that there can be 
more than two minimizing curves starting at a shock point x; their velocities are between 
v(x + ) and v(x^). This happens, e.g., at shock mergers. Since \v(x)\ = p(y' J ,y'j_ 1 ), we 
have: max x \v(x)\ < 1/2, max a . |uoo(x)| < 1/2 + max^ |y(x)|. We have seen above that the 
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points of continuity of u^, and hence of v, are exactly the points of uniqueness of minimizing 
curves. Since any limiting point of a sequence B^uq(x) — g(x) gives a velocity of a minimizing 
curve starting at (x, J), we have lim rwoo (i3™uo(x) — g(x)) = v(x), or B^uq(x) — ► Uoo(x) as 
n — > oo for any x outside of S. Suppose now that there exists another global minimizer 
{yj}- The same argument as above shows that necessarily lim^oo j/j = lim ; 7_ + _ 00 y~j = x c . 
If y~i 7^ x c for some i, then one can construct a sequence with smaller action by taking 
jjj = x c , \j\ < jo, Vj = Vj, \j\ > jo for sufficiently large jq. Such construction contradicts to 
global minimality of {yj} and proves uniqueness of the global minimizer. Statements 1-4 
are thus proved. 

Uniqueness of the main shock 

Till now we have not used the hyperbolicity properties of the fixed point P. We have seen 
above that if v is a velocity of a minimizer which starts at (x, J), then the point (x, v) belongs 
to I?( u '. Denote by s the natural parameter of length along TW, i.e. 

r( u ) = (x(s),v(s)), seR, (dx/ds) 2 + (dv/dsf = 1, (x(0),v(0)) = (x c ,0). (A.6) 

The orientation of s is fixed by the condition x(s) | x(0) as s j 0. Define (xj(s),yj(s)) = 
T- 7 (x(s),f(s)), j < and the C 2 function 



Ms) = E 

j<o 



^H> - G(y 3 



(A.7) 



The series above converges since G(x c ) = 0. It is easy to see that a point (x(s),v(s)) 
corresponds to a minimizer if and only if A(s) = ming^s), the minimum being over all s 
such that x(s) = x(s). Denote by A{x) = min^ A(s), the minimum being now taken over all s 
such that x(s) = x. Notice that for any e > there exists 5(e) > such that A(s) > 5(e) for 
all |s| > e. Fix e small enough, so that x(s) is a monotone function for |s| < e. Now, choose 
eo so small that max| s |< eo A(s) < 5(e). Then for all \s\ < eo, a point (x(s),v(s)) corresponds 
to the unique minimizer at x(s). Hence, we have shown that there are no shocks inside some 
neighborhood of x c . 

We now construct the main shock. Fix an arbitrary time t. Consider the situation on the 
space-time cylinder M. All minimizing curves approach the global minimizing curve 7 either 
from the right or from the left (see figure pi). Denote by B r (t) and Bi(t) the sets of points x 
on the circle such that there exists a minimizing curve starting at (x, t) which approaches 7 
from the right and from the left, respectively. Since minimizing curves do not intersect, B r (t) 
and Bi(t) are closed intervals, B r (t)\J Bi(t) = S 1 = [0, 1[ and B r (t) f] Bi(t) consists of just 
two points. One of them is x c . Denote the other one x ms h(t). It follows immediately from the 
construction that x ms h(t) is a shock point; moreover it is the main shock. To prove uniqueness 
consider any other shock at time t = J at point x. Then, (x,i>(x_)) and (x,t>(x + )) either 
belong both to the negative-s part of r( u ) or both to its positive-s part; that is, there exists 
s\,S2 such that (x,t>(x_)) = (x(s\),v(si)) and (x,t>(x + )) = (x(s2),v(s2)) with S1S2 > 0. 
For j sufficiently negative, both (xj(si),yj(si)) = T^(x(si),v(si)) and (xj(s2),Vj(s2)) = 
T 3 (x(s2),v(s2)) belong to an eo-neighborhood of (x(0),w(0)) where there is no shock. Hence 
a prehistory of an original shock is not longer than \j\. Statement 5 is thus proved. 

Exponential convergence to the unique solution 

Consider a small neighborhood U of the point P. It is well known (Hartman 1960; Belitskii 
1973) that if U is small enough then, inside U, T is C 1 -smoothly conjugate to a linear trans- 
formation. This means that there exists a local C 1 -smooth change of variables such that in 
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the new coordinates (X, V) the map T becomes T\ : (X, V) i— > (AA, A" 1 ^) and P is the 
origin in the coordinates X, V. Denote u n (x) = B™uq and recall that J uo(x)dx = 0. For 
arbitrary x consider a point (x, v) = (x, u n (x)) on the cylinder C and its backward trajectory 
( x (— J)> v (— j)) = T~i(x,v), < j < n. It is easy to see that there exists a neighborhood 
U\ C U such that if two points of this backward trajectory belong to U±, then all points 
in between belong to U. Notice, that there exists rii which depends only on U\ such that, 
uniformly in n, at most n\ points of the backward trajectory (x(—j),v(—j)), < j < n, are 
outside of U±. This implies that for some j±, < j± < m, and j 2 , n — n\ < j 2 < n, we have 
(x{-ji),v(-ji)) G U 1 ,(x(-j 2 ),v(-j 2 )) G U 1 . Hence, (x(-j),v(-j)) G U for all ji < j < j 2 
and j'2 — ji > n — n\. Then, there exists a constant C\ > such that the distance between 
(x(—ji),v(—ji)) and a piece of local unstable manifold inside U is less than C\X~ n . Denote 
by ]s'" , s^ 2 '[, s' 1 ^ < 0, s^ 2 ' > an interval of the parameter which corresponds to this piece of 
local unstable manifold. Then the Euclidian distance dist ((x(— j±), v(— ji)), (x(s n ),v(s n ))) < 
C\X~ n for some s n G [s^^s^ 2 ']. Denote by s' 1 ^ and s' 2 ' the values of the parameter s cor- 
responding to T ni (x(s^ 1 >),v(s^ 1 ')) and T ni (x(s^ 2 '),v (s' 2 ')), respectively. Then, there exists 
a constant C > such that dist ((x,v), (x(s n ),v(s n ))) < CX~ n for some s n G [s^ l \ s^ 2 ']. 
It might happen that s n does not correspond to a minimizer. However, A(s n ) — ► A(a;) as 
n — > cxd. More precisely, one can show that there exists a constant C 2 > such that 

A(a n ) - A(x(s n )) < C 2 nX- n , \A(x(s n )) - A(x)\ < C 2 X~ n . (A.8) 

Suppose now that Uoo is differentiable at x (recall that it is differentiable almost everywhere). 
Denote by s x the value of parameter s corresponding to the unique minimizer at x. Then, 
dx(s)/ds \s= Sx 7^ 0- Then, there exist e, 5 > such that \dx(s)/ds\ > 5, \dy(s)/dx\ < S^ 1 for 
all s e]s x -e,s x + e[. Denote by v = mm se[ ^ 1)sx _ e] y [aB+e)g (2)] dist {(x, A(x)), (x(s), A(s))} > 
0. Clearly, there exists N which depends only on v such that dist {(x, A(x)), (x(s n ), A(s n ))} < 
v for all n > N . Hence, s n £]s x — e, s x + e[ for all n > N. This and the estimates for the 
derivatives immediately imply that there exists a constant C(x) > such that \B™uq(x) — 

«oo(aOI < CWA"". 

To prove an estimate from below, notice that, for any backward trajectory 
( x (— J)) v (— j))-, < j < n, the last point (x(— n), v(— n)) cannot be too close to P if uo(x c ) / 
0. Indeed, the initial potential ipo has non-zero slope at x c since we assumed uq(x c ) ^ 0. 
Hence, it is possible to make the action smaller by moving further from x c . It is easy to show 
that there exist e(uo) > such that dist {(x(— n),v(— n)), (x c ,0)} > e(uo). It follows that 
a point (x,v) cannot be too close to (x(s),v(s)), s G [s^ 1 ', s( 2 >], i.e. there exists a constant 
c(uq) > such that |-B"«o(x) — Uoq(x)\ > c( , uo)A -n . Statement 6 is thus proved. 

Analyticity and finiteness of the number of shocks 

To prove Statement 7 we first notice that T^ n > is analytic, provided G is analytic (see Moser 
1956). More precisely, Moser 's result implies local analyticity, from which we can infer an- 
alyticity on any closed interval of r( u ) which does not contain the fixed point. Denote by 
{xi} the (at most countable) set of shock points other than the main shock at an integer 
time J. As we have already seen above, for every x» there exists an open set ]s±(i) , s 2 (i)[ of 
the parameter s, where the parameter values si(i), s 2 (i) correspond to the points (xj, v(xi_)) 
and (xi,v(xi + )), respectively. Since x« is not the main shock, si(i)s 2 (i) > 0. The non- 
intersecting property of minimizing curves implies that different intervals ]s\(i), s 2 (i)[ do not 
intersect. Also, ]si(i),s 2 (i)[c [s2msh> simsh]> where si ms h and S2msh are the parameter val- 
ues corresponding to (x ms h,v (x ms h_)) and (x ms h, , w(x ms h + )), respectively. It is easy to see 
that all s G [s2msin si ms h] — Ui(si(i), s 2 (i)) correspond to minimizers. Clearly, for all i there 
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exists s(i) G (si(i),S2(i)) such that dx(s)/ds \ s=s u-\ = 0. Suppose that there be infinitely 
many shocks. Then there exists an accumulation point Sqo for the sequence {s(i)}. It follows 
that all derivatives of x(s) vanish at Sqq. The analyticity of x(s) then implies that x(s) is a 
constant function. This contradiction proves that the number of shocks is finite. The same 
argument works for all times t, since for all t the Lagrangian map transforms any piece of 
finite length of r^ u ' into an analytic curve. Denote by T\ u = (xt(s), vt(s)) for < t < 1 the 
image of T^ u ' under the Lagrangian map at time t, where s is a natural parameter along T t n . 
Suppose that the number of non-main shocks at time t is K, so that the total number of 
shocks is K + 1. As above, denote by ]si(i), S2(i)[ the intervals of the parameter s generated 
by the i-th. shock, and by S2msh and si ms h the values of the parameter corresponding to the 
main shock. The intervals ]s\(i), S2(i)[, 1 < i < K, divide [s2msh> s imsh] into K + 1 closed 
intervals I{, 1 < i < K + 1. It is easy to see that each of those intervals corresponds to an 
analytic piece of Uoo(x,t), i.e. Uoo(x(s),t) = vt(s), s £ h. Finally, we show that the number 
of preshock events between time and 1 is finite. Suppose the number of preshocks were 
infinite. Denote by U,Si the time of the i-th preshock and the corresponding value of the 
parameter s. Denote by (t*,s*) an arbitrary point of accumulation for a sequence (ij,Sj). It 
is easy to see that all derivatives of xt*(s) vanish at point s*. This implies that xt*(s) is a 
constant function. Again, we get a contradiction, which finishes the proof of Statement 7. 
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